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ABSTRACT 

Optical and infrared interferometers definitively established that the photo- 
metric standard Vega (= a Lyrae) is a rapidly rotating star viewed nearly pole-on. 
Recent independent spectroscopic analyses could not reconcile the inferred incli- 
nation angle with the observed line profiles, preferring a larger inclination. In 
order to resolve this controversy, we observed Vega using the six-beam Michi- 
gan Infrared Combiner (MIRC6) on the Center for High Angular Resolution 
Astronomy (CHARA) Array. With our greater angular resolution and dense 
(u,v)-coverage, we find Vega is rotating less rapidly and with a smaller gravity 
darkening coefficient than previous interferometric results. Our models are com- 
patible with low photospheric macroturbulence and also consistent with the pos- 
sible rotational period of ~0.71 days recently reported based on magnetic field 
observations. Our updated evolutionary analysis explicitly incorporates rapid 
rotation, finding Vega to have a mass of 2.15lo!i5 M Q and an age 700+ 7 5 Myrs, 
substantially older than previous estimates with errors dominated by lingering 
metallicity uncertainties (Z = 0.006to!oo2)- 

1 monnier@umich.edu: University of Michigan (Astronomy), 500 Church St, Ann Arbor, MI 48109 

2 The Pennsylvania State University, University Park, PA 16802, USA 

3 Geneva Observatory, University of Geneva, Maillettes 51, 1290, Sauverny, Switzerland 

4 Sydney Institute for Astronomy (SIfA), School of Physics, University of Sydney, NSW 2006, Australia 

5 Embry- Riddle Aeronautical Univ, Daytona Beach, FL 32114 USA 

6 Centre de Recherche Astrophysique, Ecole Normale Superieure de Lyon, F-69384 Lyon cedex 07, France 
7 The CHARA Array of Georgia State University, Mt. Wilson, CA, 91023, USA 
8 The Scottish Association for Marine Science, Dunstaffnage, Oban, Argyll PA37 1QA UK 
9 National Optical Astronomy Observatory, 950 N. Cherry Ave, Tucson AZ 85719. USA 
10 Department of Physics and Astronomy, University of St. Andrews, Scotland UK 



- 2 - 



Subject headings: stars: individual (Vega), techniques: interferometric, infrared: 
stars, stars: rotation 



Introduction 



The nearby hot star Vega (spectral type AO) has been used as a photometric standard 
for millennia. While Vega's relatively narrow spectral lines (vsini ~22 km/s) suggest slow 
rotation, interferomet ric observations instead have established Veg a to be a rapid rotator 



viewed near pole-on (IPeterson et al 



of earlier spectroscopists (IGray 



1985 



2006 



1988 



Aufdenb erg et al. confirming suspicions 



Gulliver et al 



19941 ) . Rapid rotation should 



keep the surface material and stellar envelope well-mixed, leading to the conclusion that 
the observed sub-solar photospheric abundance represents the bulk composition. Lower 
metallicity l ed to a revised lower m ass estimate of ~2.15M for Vega and increased age 
-500 Myrs (TYoon et alJI^Ooi . boioh . 



Based purely on spectroscopic analysis, Takeda and collaborators (ITakeda et al.ll2008al : 
Takeda et al. I l2008bh agree that Vega is rapidly rotating but with a preferred set of parameters 
at odds with the first-generation of interferometry results. The parameters most discrepant 
are the rotational period and inclination angle, key values for modeling line profiles and un- 
derstanding its evolutionary state. Yoon et al. (2008, 2010) made the case that non-standard 
macroturbulence broadening of ~10 km/s could accommodate both the observed line profile 
shapes and the interferome try results (Take da et al. adopted 2 km/s microturbulence with 
no additional broadening). iHill et al.l ( 120101 ) carried out a similar analysis and found inter- 
mediate results for best-fitting macroturbulence and inclination angle. We refer here to this 
tension between models as Vega's inclination controversy, although one might alternatively 
refer to it macroturbulence controversy. 

A third observing method has r ecently shed new light on this touchstone system. From 
analysis of circularl y-polarized l i ght. iLignieres et al.l ( 120091 ) found evidence for a weak mag- 
netic field in Vega. iPetit et al.l (120101 ) carried out Zeeman Doppler imaging, finding a de- 
tectable weak polar field concentration (~0.6 Gauss). Although the periodic signal is indeed 
weak, there is grow ing confidence after many years of observations that a persistent signal 
at 0.71±0.03 days (IPetit et al.ll2010l : lAlina et al.l 120121 ) represents the rotational period of 
Vega. This period is substantially longer than expected from interferometry-based models 
(P~ 0.5-0.6 days) but compatible with the period range predicted using line profiles alone 
(P~ 0.7-0.9 days; Takeda et al. with no excess macroturbulence). 



In this paper, we present extensive new interferometer observations of Vega using the 
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Michigan Infrared Combiner (MIRC) on the Center for High Angular Resolution Astronomy 
(CHARA) Array. Our data have higher angular resolution than previous work and sub- 
stantially improved Fourier coverage, allowing a robust estimate of internal and systematic 



error s . Building on our recent imaging and modeling of other rapid rotators (IMonnier et al. 



20071 ; IZhao et al.l 120091 ; IChe et al.l 120111 ) , we explore a wider range of gravity darkening pre- 



scriptions. In short, we bring a new independent and critical look at the constraints interfer- 
ometers can bring to Vega, particularly cognizant of parameter degeneracies and calibration 
systematics. 



2. Observations and Data Reduction 



We have used CHARA array in conjunction with the MIRC combiner to measure visi- 
bilities (V 2 ) and closure phases (CP) of Vega across the near-infrared H-band. The CHARA 
Array was built and is operated by Georgia State University on Mt. Wilson, California. 
CHARA is the longest baseline optical/infrared i nterferometer in the world w ith six fixed 
1-m telescopes and a maximum baseline of 330m (jten Brummelaar et al.l 120051 ). 



The MIRC image-plane combiner was used for all observations presented here. Before 
2011, MIRC was used to combine four telescope beams, allowing 6 V 2 and 4 CP measure- 
ments at a time. Following a major upgrade in 2011, MIRC now combines all six CHARA 
telescopes, resulting in up to 15 V 2 and 20 CP measurements simultaneously. MIRC splits 
the H band light (Aq = 1.65/im) into eight spectral channels ~42), with absolute wave- 



l ength precision of ±0.25% based on measures of i Peg using the orbit of Konacki et al. 



2004 


.12006. 


2010; 


Che et al. 


2010, 


2012) 



Using Fourier transform techniques, the V 2 are measured, averaged and corrected for 
biases. The bispectru m is formed us ing the phases and amplitudes of three baselines that 
form a closed triangle ( iMonnierl 120071 ). Amplitude calibration was performed using realtime 
flux estimates derived from choppers (before 2010) or t hrough use of a beamsplitter following 
spatial filtering for improved performance (after 2010; IChe et al.ll2010l ). Lastly, observations 
of reference calibrators throughout the night allowed for correction of time-variable factors 
such as atmospheric coherence time, vibrations, differential dispersion, and birefringence in 



2007; 


Zhao et al. 


2009; 


Che et al. 


2011) 



For this work we have evolved our calibration model to better account for systematics. 
Firstly, we include two types of calibration error for V 2 - multiplicative errors associated 
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with the transfer function and additive errors associated with correcting biases at low fringe 
or bispectrum amplitude. Based on calibrator studies, the former has been estimated to be 
20% (6.6%) for 2007 (2012) data while the additive systematic error is AV 2 =2 x 10~ 4 (for 
both epochs). For triple amplitudes ("T3amp"), the corresponding multiplicativ e errors are 



30% ( 10%) for 2007 (2012) and additive errors are 1 x 10~ 5 . A detailed study bv lZhao et al. 



( 120111 ) suggests CPs have an error floor of 1° for the observing modes adopted here. To avoid 
our model fits being trapped by systematics, we also include two new types of CP errors 
associated with low signal-to-noise ratio (SNR) data near visibility null crossings. Because 
correlated camera readout noise dominates the CP measurements at low SNR, we enforce 
minimum CP errors when the SNRT3amp~l- I n addition, we account for finite time-averaging 
and spectral bandpass effects by including an error term proportional to ACPa across each 
spectral channel. Formally, these two terms are only important right at the null crossings 
and appear in the following noise floor formula: <x CP > MAX(™^ ^,0.2 ACP A ). 

(oiN -tlT3amp/ 

Here we present data for Vega from 3 nights in 2007 (MIRC4) and 2 nights in 2012 
(MIRC6). Data on four additional nights of MIRC4 were recorded in 2007 and 2010 but 
were discarded due to calibration problems. Table [1] includes detailed observing information 
includin g the calibrators and their adopted sizes; reduced data are available in OI-FITS 



format (IPauls et al.l 120051 ) upon request 



Inspection of CPs show nearly all values to be at zero or 180° as expected for a point- 
symmetric intensity distribution. Figure 1 shows the visibility data and the (u,v)-coverage 
(inset) for our datasets, split into three chunks of similar quantity: 2007, 2012 Jun 9 and 
2012 Jun 13. The data was azimuthally-averaged and compared to uniform disk and power- 
law limb-darkened disk ( / = I fi a ) models. As expected, the data are not consistent 
with a uniform disk and we find a best-fit limb-darkened diameter of 3.324 milliarcseconds 
(mas) with power law a ~ 0.227 - more limb-darkened than expected for a non-rotating 
star (a ~0.11, Kurucz). We note variation between epochs due to calibration errors. Mea- 
suring precise limb-darkening requires controlling systematics at the few % level, a goal for 
CHARA/MIRC but one that still proves challenging to attain during all observing condi- 
tions. 

Internal diagnostics demonstrate that the best calibration is from 2012 Jun 13. This 
night had the best flux calibration, the greatest on-source integration time, and also employed 
the maximum number of simultaneous telescopes. The 2007 data, while extensive, were taken 
before our photometric channel upgrade and suffer from larger calibration errors. For the 
detailed modeling in the next section, we limited our fits to the 2012 Jun 13 dataset. 
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3. Modelling 



Our group has been a leader in the field of modeling rapid rotators, based largely on 
our unique and extensive interferometry data from CH ARA/MIRC. Our series of papers 
( IMonnier et al.l 120071 ; IZhao et al.l 120091 ; I Che et al.l 1201 f I ) contains the first images of main- 
sequence stars beyond the Sun and has determined precise stellar parameters for rapidly 
rotating stars from early F to late B. The basic physical model consists of a star with 
uniformly-rotating surface layers, distorted by centrifugal forces acting under point-gravity 
with a surface temperature following the gravity-darkening law T = T^oi^T- 2 — where g 



is the effective surface gravity. Based on our full dataset, IChe et al.l ( 120111 ) argued that 
the observed grav ity darkening deviates from the canonical value of (3 = 0.25 advocated by 
von Zeipell (Il924al lbl) and instead we find empirically a lower characteristic value of /3 = 0.19. 
For our work here, we will again consider a range of possible coefficients. Details of our cal- 
culati ons can be found in these earlier papers which followed the method of lAufdenberg et al. 
( 120061 ); the full list of independent model parameters is shown at the top of Table [2J 



In ord er to car r y out the full calculation, we had to assume a few physical parameters. 
We used a iKuruczl (119791 ) grid (see |http: / /kurucz.harva rd.edu/grids.html ) for a [Fe/H|=- 
0.5 ( sub solar) plane-parallel atmosphere (recent metalli city determinations by lYoon et al. 
2010). We also used the Hipparcos distance o f 7.68 pc (jvan Leeuwenl 120071 ) and fixed the 
model mass to 2.15 M as recommended from Yoon et al. ( 201ol . and consistent with final 
results presented here). 

X 2 -minimization was used to const rain our model paramet e rs by e mploying the ensemble 



Markov Chain method descr ibed bv Foreman- Mackev et al.l ( 120121 ). based on the affine- 
invariant strategy outlined in Goodman fc Weearel ( 2010 ). We used 1000 walkers seeded by 
uncorrelated random distributions of the entire relevant parameter space. The distribution 
reached convergence typically within 25 steps although we calculated 75 steps before freezing 
the distribution for error analysis. The statistical weight of the V 2 , T3amp and CP fits were 
reduced by the number of spectral channels (= 8) because of strong internal correlations. 
In addition, we down-weighted the contribution from the V 2 and T3amp by an additional 
factor of two, since these quantities are not independent of eac h other (i.e., T3amps ca n be 
derived from the V 2 ). The observed V-magnitude (0.03 mag, iMermilliod et al.l Il997l ) and 
H-magnitude (0.00 mag. lKidger &; Martm-Luisll2003l ) photometry (adopting 5% photometric 
errors to account for both measurement and zero-point uncertainties) was incorporated as 
a statistical prior during the Markov Chain calculation. Our final error bars combined 
the errors from the Markov Chain in quadrature with errors due to calibrator diameter 
uncertainty - this was done by repeating the above Markov Chain for different assumed 
calibrator diameters drawn from the expected range of values. These two error sources were 
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often similar in magnitude, although some parameters have uncertainties dominated by one 
or the other. 

Table [2] shows the full results of three separate parameter studies, labeled Model 1, 2, 
& 3. Becaus e of the severe d egeneracy between uj and (3 for pole-on rapid rotators (see 
discussion in IZhao et al.ll2009l ). we chose to fix the gravity darkening coefficient for Model 1 
to be = 0.25, consistent with the classical von Zeipel value allowing for compa rison with 
previo us work. For Model 2, we fixed (3 = 0.19 equal to the recommended value of IChe et al. 
( 120 111 ). Before describing Model 3, we first discuss the results of Models 1 & 2. 



Most of the model parameters are consistent (within uncertainties) between Models 
1 & 2. Model 1 prefers a slower rotational rate (longer periods) than Model 2, which 
was expected since a faster rotation state is needed to compensate for the weaker gravity 
darkening to maintain the center-to-limb intensity profile. Note that the x 2 f° r Models 1 
and 2 are practically identical (xl = 0-89 and 0.88, respectively), showing the near perfect 
degeneracy in ft-u space for a pole-on system. While our models do prefer low inclination 
angles (i ~ 4.5°), in agreement with ea rlier interferometric models ([Peterson et al.l 12006 



Aufdenberg et al.ll2006t lYoon et al.ll2010l ). the determination of our uncertainty is nearly 5x 
larger than those derived by earlier workers. This is true for other parameters too - our 
analysis takes a much more modest view of parameter uncertainties which largely eliminates 
the strong conflict with spectroscopic analyses. For instance, we find v sinz = 15 ± 4 
km s _1 (no macroturbulence) which is more compatible with the observed 22±2 km s _1 . The 
most notable statistically-significant disagreement with previous determinati ons is uj / uj cr \t' 
for {3 =0.25, we fi nd uj / cu C r\t= 0.77 ± 0.0 5, not compatible with 0.93±0.02 (IPeterson et al. 



20061), 0.91 ±0.03 (lAufdenberg et al.ll2006l ). or 0.876±0.006 (lYoon et al.ll2010h . We speculate 
that the more limited datasets of these workers led to underestimates of systematic errors 
although we can not rule out certain physical explanations (time variability, spots, non- 
standard gravity darkening). 

Constrained by only interferometry data, our parameters of Models 1 & 2 span a larger 
range than earlier estimates, limiting our constraints on key stellar prop erties. We can 



reduce our errors by including constraints on the period (0.71±0.03 days; iPetit et al.l 12010 
Alina et al.l 120121 ) and vsini (22±2 km s _1 , Takeda et al. 2008ab) as statistical priors 



during the Markov Chain calculation. We call Model 3 our "Concordance Model," a set of 
parameters that agrees with CHARA/MIRC interferometry, the SED, vsini, and rotational 
period estimates. Table 2 contains these results, showing much smaller error bars with only a 
minimal increase in normalized x 2 (0.89 to 0.90). In this model we allowed gravity darkening 
to be free and find (3 = 0.231±0.028, a bit higher than the 0.19 pr eferred bylche et al.l (boilh 
but interestingly consistent with the re-parametrization of (3 by lEspinosa Lara &: Rieutord 
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( 120 111 ), who argue that /3 depends on the rotation rate, matching 0.25 only for slow rotators 



becoming smaller as the star spins up. 

As a side note, we found that the basic parameters of the Concordance Model can 
be determined simply and robustly without a complicated calculation. First, the projected 
equatorial diameter can be deduced from basic visibility fitting (see Figure 1) to be 3.32 mas. 
This diameter can be turned into an equatorial velocity (t> eq = 195 km s" 1 ) using the 
estimated period and distance. And finally, the inclination angle must be ~6.5° to match 
the observed vsmi. Indeed, these "back-of-the-envelope" estimates match quite closely the 
Concordance Model results in Table 2. 

One way to vie w these result s is to compare our parameters with the family of SED 



solutions outlined by lTakeda et al.l ( l2008af ). Figure 2 shows our interferometric-based models 
(along with those of Peterson et al. and Aufdenberg e t al. ) plotted in cj-inclination space 
next to Takeda's Models #1-93 • Takeda et al. ( 2008a ) went further and used the line pro- 



file shapes to select best Models #3, 4, 5, corresponding to ~7° inclination. We see our 
Models 1 & 2 are just consistent with the Takeda results at the ~l-sigma level while our 
Concordance Model strongly selects Takeda Model #5 as the optimal choice. Indeed, this 
diagram further reinforces that a true concordance does exist between the SED and line pro- 
file fitting of Takeda, the putative rotational period from Petit et al. and the CHARA/MIRC 
interferometric observations in the near-infrared. 

We end this section with a note of caution. The deviation from centro-symmetry on 
the surface of Vega is quite subtle, amounting to a pole offset of just ~ 0.2 milli-arcsecond, 
roughly 5x smaller than the fringe spacing from our longest baseline. Model-fitting different 
MIRC epochs can yield pole PAs as diffe rent as 90°, and result s from our best epoch (see 



Table 2) are discrepant with results from iPeterson et al.l ( 120061 ). We have searched exten 



sively for the explanation for the fragile constraints on the pole PA, inclu ding physical causes 



(faint close-by companion, magnetic spots, non-radial pulsations; see also lRogers et al.ll2012f ) 
and calibration-related problems (fringe cross talk, detector noise, bispectrum bias, time- 
averaging). After an exhaustive series of tests, none of these hypotheses could convincingly 
explain the variations. We urge follow-up observations, especially at visible wavelengths 
where gravity darkening effects are strongest. Fortunately, the conclusions from our work 
here depend mostly on V 2 and not CPs, and the pole PA is not of paramount physical 
importance. 



*Takeda's curve was slightly shifted here to account for fact that we use 2.15 M instead of 2.3 



Sini ncw = J A '^ da Sinixakoda 
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Discussion and Conclusions 



We can now use our modeling results to assess the evolutionary state of Vega. We 
have used the most recent evolutionary tracks from the Geneva grou p that explicitly incor - 
porates the effects of rapid rotation (Georgy et al. 2013, submitted; lEkstrom et al.ll2012l ). 
We considered metallicity range Z = 0.006^q;qq2 appropriate for \Fe/H] = —0.5 under the 



range of currently conside red che mical abundances of the Sun (lAnders fc Grevessd 11989 



Asplund et al.ll2005l 120091 ) - note lYoon et al.l (120101 ) recommend Z = 0.009 corresponding 



to the upper range we considered. Because rapid rotation makes a star's posi tion on the 



tradi tional H-R diagram (L vs. T eff ) to be viewing angle dependent (see e.g., IZhao et al. 



20091 ). we instead present stellar evolutionary tracks in units of total bolometric luminos- 
ity and stellar polar radius. Figure 3 shows our modified Hertzsprung-Russell diagram for 
Z = 0.006 including the effect of rotation. Our best models have to / u; cr i t ~0.8 and we 
show these isochrones in the figure. We conclude that Vega has a mass of 2.15^^5 M and 
age of 700+j| million years for Z = O.OOGj^ooi) with errors dominated by the metallicity 



assum ption not random errors. While our mass estimate is similar to those of lYoon et al. 



(120101 ) . our age estimate is significantly higher due mostly to including the effect of rotation 
and less so because of the lower mean metallicity Z we have adopted. 

In conclusion, we have presented modeling of the photometric standard star Vega using 
new interferometric data from CHARA/MIRC. The large quantity and high angular reso- 
lution of our data allow for precise constraints on the geometry and surface temperatures 
of Vega. We find Vega rotating more slowly th an previous i nterfe rometer results, consistent 
with the putative rotation period observed by lAlina et al.l ( 120121 ) and compatible with the 
observed line profiles without excess macroturbulence. The differences with previous inter- 
ferometry results could be from under-estimates of errors in earlier work or may suggest 
subtle deficiencies in the physical models. Our "Concordance Model" and its placement on 
a new H-R diagram represent the best global model for Vega to date but there is still room 
for improvement. In addition to confirmation of the rotation period through photometry, we 
recommend additional visible-light interferometry data spanning the first 3 visibility lobes 
with <5% precision on V 2 to definitively establish the tilt angle of the pole and to pinpoint 
the true level of gravity darkening. 



JDM thanks Deane Peterson, Yoichi Takeda, and Pascal Petit for sharing their insights 
into Vega. The CHARA Array is currently funded by the National Science Foundation 
through AST-1211929 and by the Georgia State University. Funding for the MIRC combiner 
came from the University of Michigan and observations were supported through NSF grants 
AST-0352723, AST-0707927, and AST-1 108963. STR acknowledges partial support from 
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NASA grant NNH09AK731. This research has made use of the SIMBAD database, operated 
at CDS, Strasbourg, France, and NASA's Astrophysics Data System (ADS) Bibliographic 
Services. 
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Table 1. CHARA/MIRC log for Vega observations 



Date 


Interferometer 


Number of 


Number of 


Calibrator 


(UT) 


( Configuration) 


V 2 


Closure Phases 


Information 


2007 Jul 5 


S1-E1-W1-W2 


168 


104 


a Cyg a , T Peg b 


2007 Jul 8 


S1-E1-W1-W2 


96 


64 


7 Lyr c T Peg 


2007 Jul 13 


S1-E1-W1-W2 


144 


96 


f Cyg 


2012 Jun 9 


S2-S1-E1-E2-W2 


200 


144 


HD 167304 d 


2012 Jun 13 


W1-S2-S1-E1-E2-W2 


560 


640 


7 Lyr 



a Adopted a Cyg UD diameter 0.54 ± 0.02 mas llBarnes et al.lll978h 

b Adopted T Peg UD diameter 0.99 ± 0.02 mas (new CHARA/MIRC measurement) 

c Adopted 7 Lyr UD diameter 0.737 ± 0.015 mas based on independent measurements by 
CHARA/MIRC (UD H = 0.723 ± 0.025 mas) and CHARA/PAVO (UD H = 0.744 ± 0.019 mas 
derived from UD LDD = 0.755 ± 0.019 mas) 

d Adopted HD 167304 UD diameter 0.69 ± 0.05 mas llBonneau et alj|2006h 
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Table 2. Modeling Results for Vega 



Model Parameters' 1 


Model 1 


Model 2 


Model 3 




= 0.25 


/3 = 0.19 


Concordance 15 


Inclination (°) 


4.5 ± 1.3 


4.0 ±1.5 


6.2 ±0.4 


Pole position angle (East of North ) 


-57 ±7 


-57 ±6 


-58 ± 6 


Tpoi(K) 


10120 ± 140 


10130 ± 140 


10070 ± 90 


Rpoi(Re) 


2.42 ±0.05 


2.31 ±0.06 


2.418 ±0.012 


Rpol(mas) 


1.47 ±0.03 


1.40 ±0.04 


1.465 ±0.007 


w / w crit 


0.77 ±0.05 


0.86 ±0.04 


0.774 ±0.012 


/3 


0.25 (FIXED) 


0.19 (FIXED) 


0.231 ±0.028 



Derived Parameters 



Tcq(K) 


8870 ± 200 


8740 ± 190 


8910 ± 130 


Req( Rq) 


2.726 ±0.007 


2.728 ± 0.007 


2.726 ±0.006 


Req(mas) 


1.651 ±0.004 


1.652 ±0.004 


1.651 ±0.004 


v / v crit 


0.58 ±0.05 


0.68 ±0.05 


0.581 ±0.012 


Bolometric luminosity I/bol(L©) 


47.1 ±2.7 


44.8 ±2.6 


47.2 ±2.0 


Apparent luminosity L app (L©) 


58.8 ±2.7 


58.8 ±2.7 


58.4 ±2.2 


Apparent effective temperature T^ p ( K) 


9680 ±110 


9670 ±110 


9660 ± 90 


Surface-averaged T e g ( K) 


9350 ± 110 


9310 ± 110 


9360 ± 90 


usini ( km s —1 ) 


15.1 ±3.6 


16.2 ± 4.8 


21.3±1.2 b 


Rotational period (days) 


0.71 ±0.07 


0.59 ±0.06 


0.71 ±0.02 b 


Model V magnitude 3. 


0.032 ± 0.028 


0.032 ± 0.028 


0.035 ±0.023 


Model H magnitude a 


0.029 ±0.013 


0.030 ±0.013 


0.031 ±0.011 


Model mass (M Q ) C 
Model age (Myrs) c 


Joint Result: 2.15^55 "5 
Joint Result: 7007j| 


Summary of x 2 Results 


Total xl (AWa = 1842) 


0.89 


0.88 


0.90 


Vis 2 xl (AWa = 560) 


1.45 


1.45 


1.44 


T3amp xl (N DATA = 640) 


0.57 


0.57 


0.58 


CP xl (AWa = 640) 


0.83 


0.81 


0.84 



a O ther parameters: dis tance 7.68 pc jvan Leeuwcn 2007), V mag 0.0 3±0.05 llMerm illiod et alj 
Il997l) . H mag 0.00±0.05 llKideer fc Martm-Luisll2003l) , \Fe/H] = -0.5 jYoon et al.ll2008l) . 



Th e Concordance Model incorporated t he observed vsin i = 22±2 km s'MTakeda et alj 
2008a) and period estimate 0.71±0.03 days iPetit et al.ll2010l ; lAhna et alj|2012l') into the fit as 
priors. 

c Geneva stellar evolutionary tracks (Georgy et al. 2013, submitted: lEkstrom et al1|2012| ) were 
used assuming ui / u) C rit= 0.8 and covering the range Z = 0.006^q'qq2- 



- 14 - 




Fig. 1. — This figure shows three epochs of visibility observations using CHARA/MIRC. The data are not 
consistent with a uniform disk fit (top line, 3.26 mas). The power-law limb-darkened disk fit (bottom line, 
3.32 mas, a = 0.23) shows twice the level of expected limb-darkening, suggesting strong gravity darkening. 
The (u,v) coverage for each epoch is shown in each panel as an inset box (±350 meters). 
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Fig. 2. — This figure shows the oj / u> CI it vs. inclination angle results for our three models from Table 2. 
The ellipses show the 1-sigma confidence intervals. We also incl ude the prev i ous in terferometric modelling 



(2006). 



Aufdenberg et al 



results from iPeterson et al 
solutions (Takeda Models #1-9) presented in iTakeda et al 



(2006), and Yoon et al.l (|2010l ) as well as the family of 



|2008al adjusted here for M* =2.15 M ). 
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Fig. 3. — This modified Hertzsprung- Russell (H-R) diagram plots the total bolometric luminosity vs. 
the stellar polar radius for our best-fit Models 1, 2, & 3 (squares), including the effect of the calibrator 
size uncertainty (±lcr). The stellar evolutionary tra cks are based on the most recent Geneva models that 
incorporate rotation (Georgy et al. 2013, submitted; lEkstrom et al.ll2012f ). and we show isochrones for cj J 
w cr it= 0.8 (dotted line). 



